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We study how a cosmological bounce with a Type IV singularity at the bouncing point, can 
be generated by a classical vacuum F{G) gravity. We focus our investigation on the behavior 
of the vacuum F{G) theory near the Type IV singular bouncing point and also we address the 
stability of the resulting solution, by treating the equations of motion as a dynamical system. In 
addition, we investigate how the scalar perturbations of the background metric evolve, emphasizing 
to cosmological times near the Type IV singular bouncing point. Finally, we also investigate which 
mimetic vacuum F{G) gravity can describe the singular bounce cosmology. 
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I. INTRODUCTION 

The observational data that came into play the last twenty years, have made compelling the consistent description 
with a theory that can harbor all the different observationally verified phenomena. Two are the main observations 
that need to be described theoretically, the late-time acceleration verified in the late 90’s [l|, and the early-time 
acceleration. With regards to the latterjthere is much way to be covered until we conclude whether inflation ever 
existed. However, the latest Planck data pose severe constraints on inflationary models, and indicate which features 
should a consistent theory of inflation have. Modified gravity theories offer a promising and solid theoretical framework 
that can consistently describe late-time and early-time acceleration, for reviews on this vast issue see Q . Among the 
most promising are the F{R) theories of gravity, which also offer the possibility to describe simultaneously early and 
late-time acceleration Q. The F{R) gravity is the simplest modification of Einstein-Hilbert gravity, since instead of 
having simply the Ricci scalar R in the Lagrangian, a function of R appears. In four dimensions, instead of this simple 
modification of Einstein-Hilbert gravity, there also exist other theoretical descriptions that are also promising, like 
the F{G) theories of gravity, with G being the Gauss-Bonnet invariant G = R^ — AR^j^^R^’' + Rfj.upaR^''^'^, and with 
Rpi^, Rpupa standing for the Ricci tensor and the Riemann tensor respectively. Although in principle the resulting 
equations of motion are expected to have fourth order derivatives of the metric tensor, it turns out that the theory 
contains only second order derivatives, and therefore it is rendered not too complicated to be studied. Particularly, it 
has been shown that within the context of F (G) theory, late-time acceleration can be achieved [Ml. For informative 
reviews on F(G) theory of gravity see Q. 

On the other hand, an appealing alternative to the standard inflationary paradigm, comes from the bounce cos¬ 
mology theories |14l - l^ . In this kind of theories, one of the most severe drawbacks of the inflationary paradigm, the 
initial singularity problem, is absent, and it is conceivable that this feature makes them quite appealing. The initial 
singularity is a crushing type timelike singularity, and at the point it occurs the geodesics cannot be continuously 
be extended, that is, geodesics incompleteness occurs. The singularity theorems of Hawking and Penrose [ll fully 
describe these singularities, and much work has been devoted in studying spacelike singularities. In cosmology how¬ 
ever, most of the singularities are timelike, and the most severe of these are the initial singularity and the Big Rip 
[^ . With regards to the latter, it is the most severe of the four types of finite singularities classified in [2^. Also 
sudden singularities were also studied in [27l - [^ . Among the four types of finite-time singularities, the most mild 
from a phenomenological point of view, is the Type IV singularity, at which no geodesic incompleteness occurs. For 
recent studies on this kind of singularity see |3lM36j | . The most phenomenologically interesting feature of theories with 
Type IV finite time singularities is that the Universe may smoothly pass through [ssl these timelike singularities, 
without any catastrophic consequences. In some cases, there might be observational evidence or indication of the 
passage of the Universe through a Type IV singularity [s^, but it is important that this is not catastrophic and it 
always is a smooth passage (^ . 

In view of the aforementioned interesting features of the bouncing cosmologies, in this paper we aim to study a 
specific bounce which contains a Type IV singularity at the bouncing point. In this way, the cosmological system 


v.k.oikonomoul979@gmail.com 



2 


is free from the initial singularity, but at the same time we investigate the implications of such a mild finite time 
singularity. Particularly, we shall investigate how such a cosmology can be described by a classical vacuum F(G) 
gravity, with special emphasis being given on how the F{G) gravity behaves near the Type IV singular bouncing 
point^. Also, we investigate how the resulting solution behaves, by checking the stability of the resulting equations, 
when these are viewed as a dynamical system. In this case we examine if the solution is the final attractor of the 
system. It is conceivable of course that since the Universe passed through the singular point smoothly, the resulting 
F(G) gravity should not be the final attractor of the theory, so in some sense instability of the dynamical system is 
anticipated. In addition, we shall investigate how the scalar perturbations of the background metric behave for the 
case of the resulting vacuum U(G'), near the Type IV singularity. Finally, we shall investigate which vacuum F{G) 
gravity generates the same singular bounce we are discussing, but in the context of mimetic F(G) theory, which was 
developed in [s^. Again the focus will be for cosmic times near the singular bouncing point. We believe that our 
work will provide some new information on the behavior of classical modified theories of gravity near a Type IV 
singular bouncing point, a study which combines bouncing cosmology with a mild singularity at the bounce, and also 
a classical description with an F{G) gravity. 

This paper is organized as follows: In section II we briefly review the basic features of a Type IV singular bounce 
cosmology, and in section III by using known cosmological reconstruction techniques, we investigate which vacuum 
F(G) can describe the Type IV singular bounce, by focusing on the behavior near the bouncing point. Also, we 
discuss the possible connection of the resulting F{G) theory with other viable F{G) theories. In section IV we study 
the stability of the resulting solution we found in section III, by treating the system of equations of motion as a 
dynamical system. To this end, we rewrite the system of equations of motion in terms of new variables, which make 
the study more clear from a physical point of view. In section V we investigate how the scalar perturbations of the 
background metric behave, for the vacuum F(G) theory we found in section III, emphasizing again on the behavior 
near the Type IV singular point. Finally, in section VI we investigate which mimetic F{G) gravity can describe the 
Type IV singular bounce cosmology, focusing again on the behavior near the Type IV singularity. The conclusions 
along with a discussion follow in the end of the paper. 


II. A BRIEF DESCRIPTION OF THE SINGULAR BOUNCE 


As we already mentioned, bouncing cosmology [l3 - [^ is an appealing alternative scenario to the standard in¬ 
flationary paradigm, with the most attractive feature of bouncing cosmology being the fact that there is no initial 
singularity, and there exists also the possibility of successfully describing early-time acceleration [13, HI. It has been 
recently shown however, that other types of milder singularities [^. l32l-l^ 1^ might occur during the cosmolo gica l 

evolution, without having the catastrophic consequences of the crushing type singularities, like the Big Rip one (^ . 
Particularly, in Ref. [s^ it has been demonstrated that a Type IV singularity [2^ l32l - [3^ may occur at the bounce 
point of a general bouncing cosmology, having interesting consequences. We shall briefly describe this possibility in 
this section, in order to render the presentation self-contained. 

Firstly, a cosmological bounce can be separated in two evolution eras, the contraction era, during which the scale 
factor decreases (a < 0), and the expansion era, during which the scale factor increases (d > 0). In between the two 
eras, and after the contraction era, the Universe reaches a minimal radius, where d = 0, and this is the reason that 
the bouncing scenario is free of the initial singularity. In principle, the bouncing point, that is, the point at which the 
bounce occurs, can freely be chosen, so we assume that the bounce occurs at t = ts, so when t < tg-, the Hubble rate 
F[ = d/a is negative (since d < 0), while for t > tg li is positive (since d > 0), and of course at the bounce it becomes 
equal to zero H{ts) = 0. We shall consider the bouncing cosmology with scale factor. 


i{t) = 




( 1 ) 


with a(tg) = 1, and /o an arbitrary parameter. In addition, the parameter e is assumed to be e < 1, so the bouncing 
cosmology of Eq. o is assumed to be a deformation of the well known symmetric bounce [2^, 


a{t) 




( 2 ) 


From Eq. (P) it easily follows that the Hubble rate is equal to. 


H{t)=2{l+e)fo{t-tgf^+\ 


( 3 ) 


^ From this point, when we refer to the bouncing point or the singular point, we refer to the same point. 
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For simplicity we introduce the following variables, 

/3 = 2(l + e)/o, a = 2e + 1, (4) 

so the Hubble rate becomes, 

H{t)=l3{t-tsT . ( 5 ) 

According to the classification of finite time singularities, the following types of singularities occur a.t t = tg which is 
the bouncing point, 

• q; < — 1 corresponds to the Type I singularity. 

• — 1 < q; < 0 corresponds to Type III singularity. 

• 0 < a < 1 corresponds to Type II singularity. 

• a > 1 corresponds to Type IV singularity. 

so the relevant to us case is when a > I, and since we want to render the bounce of Eq. o a deformation of the 
symmetric bounce 0, we further assume that 1 < a < 2 (or equivalently 0 < £ < ^), something which implies that 
the second derivative of the Hubble rate ([S]) diverges. As shown in [^, this can have interesting phenomenological 
consequences. Finally, we need to make sure that the Hubble rate and the scale factor never become complex, so the 
parameter a needs to be appropriately chosen, so we make the choice a = 13/11, but in general a = 2 m+i ■ ^ 

detailed account on that see [s^. 


III. SINGULAR BOUNCE FROM F(G) GRAVITY 


Having described the cosmological bounce with a Type IV singularity at the bouncing point, in this section we 
shall investigate how this type of cosmological evolution can be described in terms of a vacuum F(G) gravity. Special 
emphasis shall be given in the F(G) gravity that describes the bounce near the Type IV singularity, since the general 
problem is rather difficult to address, due to lack of analytic solutions of the corresponding differential equations. 

In order to find the F{G) gravity that describes the bounce near the singular point, we shall use some very well 
know reconstruction techniques for F{G) theories of gravity dOl. For a similar technique to the one we shall use 
here, see [l^. The Jordan frame F{G) gravity action in the absence of matter fluids (vacuum F{G)), is equal to 

53 , 


S=^ Jd^xV^iR + FiG)) , (6) 

where with Mpi = 1.22 x 10^®GeV. By varying with respect to the metric, the gravitational equations 

read, 

- lof^^FiG) - ( - 2RRp, + AR^pRP - Rpc.pRpa) E'(G) (7) 

- 2 (V^V,F'(G)) R + 2gp, (□E'(G)) i? - 4 (□F'(G)) Rp, + 4 {VpW,F'{G)) RP 
+ 4 (VpV,.F'(G)) RPp - dgp, (VpV„F'(G)) RP'^ + 4 (VpV„F'(G)) g^Pg^^Rp^.p = 0 
where the Gauss-Bonnet invariant expressed in terms of the Hubble rate, equals to, 

G = 2AH‘^ {h + . (8) 

Assuming a flat Friedmann-Robertson-Walker (FRW) background, with line element, 

ds^ = — dt^-I-a^(t) ^ dx^ , (9) 

i 

the gravitational equations of Eq. © take the following form, 

6H^ + F{G) - GF'(G) + 2dH^GF''{G) = 0 
AH + 6H^ + F{G) - GF'{G) + IQHG [h -f 77^) F"(G) 

-1- 8H‘^GF''{G) + 8H^G‘^F'''{G) = 0. 


( 10 ) 
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The reconstruction technique presented in i[ni> employs the use of an auxiliary field denoted as (p, which as was 
shown can be identified with the cosmic time t, that is p = t. By introducing two proper functions of t, namely P(t) 
and Q{t), the Jordan frame action of Eq. ([6|) becomes, 


5 = 




(11) 


and by varying it with respect to t, we obtain, 


dQ(t) ^ ^ 


dt dt 

Having Eq. (TT^ at hand, solving it with respect to t = t{G) and by substituting this in the following expression 


( 12 ) 


F{G)=P{t)G + Qit), (13) 

we finally have the F{G) gravity. It is therefore obvious that the functions P{t) and Q{t) play an important role 
in the determination of the F{G) gravity, so we now derive the differential equations that yield these functions. By 
combining Eq. m and the first of the equations appearing in Eq. (US, we obtain the following differential equation, 

d P 

Q{t) = -6H^{t)-2dH^t)—, (14) 

and finally combining Eqs. m and da, we obtain the following differential equation, 

+ 2Hit) (2Hit) - H\t)) ^ + H{t) = 0 . (15) 

When solved, Eq. (ITSl) determines the function P{t), and therefore Q{t), so upon substitution of P{t) in Eq. (fCT) . we 
get Q{t). Finally, from Eq. (IT^ we get the function t = t{G), and by substituting that in Ea. dT^ we obtain the final 
form of the F{G) gravity. We now apply this method in order to find the F{G) gravity which describes the bounce 
near the bouncing point. 


A. F{G) Gravity Near the Bouncing Point 


The general problem of finding the F{G) gravity for the singular bounce with Hubble rate (O, is rather difficult 
to deal with analytically, so we shall focus on finding the F{G) gravity near the bouncing point t = tg, which recall 
that is the point where the Type IV singularity occurs too. Consequently, at some point, we shall specify our analysis 
around the singularity. 

For the Hubble rate of Eq. ®, the differential equation (1151) that yields the function P{t) can be shown that it 
becomes. 


OR P 

— (^ - + 4(i - tsTP + 1 = 0, 

a dt^ 


which can be solved analytically to yield, 


m = - 


(t - {{t - GY - 2(i - GYa - 2YGi + 2aYGi) 

2{-l + a){-l + 2a)Y 


+ C 2 , 


(16) 


(17) 


so by substituting Eq. (IT71) in Eq. (ITT)) , we obtain the function Q{t) which appears in Appendix A, since it is too 
complicated to quote it here. By using the resulting expressions for the functions Q{t) and P{t), the final form of Eq. 
OED becomes, 


l- 2 a (^42.3aQ,^3 _ UG^P) - Gx {x°^ - 2C'i/3)) 

2^ 


(18) 


where we have set x = t — tg ior simplicity. It is obvious that the above equation is rather difficult to be solved 
analytically, so we shall simplify it by keeping the dominant terms in the limit cc —> 0, which corresponds to the limit 
t ^ tg. Therefore, by taking the limit a: —>■ 0, then Eq. m becomes, 

GiGx-^°‘ - 2iCix-^+°‘aP^ = 0 , 


(19) 








5 


which yields, 

a) = ^ , (20) 

(24q,/33)T33^ 

Then by substituting (pH)) in P{t) and Q{t), and by using Eq. (IT^ . we obtain the resulting expression for the F{G) 
gravity near the singular point, which is. 


F(G) = C 2 G + AG^^ + BG^^ , 


( 21 ) 


where the coefficients A and B are given in the Appendix A. We can give however a simpler form by exploiting the 
fact that we are interested in the limit t —)> ts- The Gauss-Bonnet invariant of Eq. (|8l), is written in terms of the 
variable x = t — tg we introduced earlier, as follows, 

G = 24a;"^+3“a/33-f 24x^“/3^, (22) 

from which it is obvious that as a: — >■ 0 (or equivalently t — >■ tg), the Gauss-Bonnet invariant tends also to zero. Hence, 
by keeping the most dominant terms from the F{G) gravity of Eq. (HI]), we get the small G limit of it, 

F{G)-G2G + BG^^ . (23) 


Therefore, in the small G limit, or equivalently near the Type IV singularity, which we chose to be the bouncing 
point, the F{G) gravity that can generate the Hubble rate dS]) is approximately described by the expression of Eq. 
(1231) . Since we are discussing for cosmological times near the bouncing point, it is worth examining how the evolution 
of perturbations behave for this F{G) model. This will be done in detail in a later section. 


B. Connection with Other Viable F{G) Gravities and Possible Late-time Behavior 

In the previous section we investigated which F{G) gravity theory can successfully describe the singular bounce 
cosmology of Eq. m, near the singular bounce. However the complexity of the resulting differential equations forced 
us to find an approximate solution, with the final F{G) being that of Eq. (1^51) . or can be further simplified neat the 
bouncing point, as 

F{G) ~ HG^Tw . (24) 

Now, it would be interesting to ask how the complete classical F{G) gravity would look like. A possible answer is that 
the complete classical F{G) gravity description would be one of the four possible forms of F{G) gravity that lead to 
finite-time singularities, first studied in [1^ and further investigated in [^,1^. Particularly, the F{G) gravities that 
lead to singularities are of the following form ElSIill, 


P(C\ _ aiG^ + bi 

(25) 

FIG) + 

a2G- + b2 ’ 

(26) 

E(G) =a3G”(53G'" + I), 

(27) 

= u,G- + h' 

(28) 


where the parameters Oi, and bi with z = I, 2, 3, are arbitrary real constants. Obviously, the F{G) gravity we found, 
which for small G is given by Eq. (El, can be the limiting case of the above F{G) gravities, which lead to finite 
singularities for some values of the parameters. For example the F{G) gravity of Eq. (EH)) in the small G limit behaves 
as ^ ^G"^, which is clearly similar to the one we obtained in Eq. El- Note that in the large G limit, the F{G) 
gravity (l28l) behaves ^ so the late-time behavior of this F{G) gravity is described by a power-law F{G) 

function and as it was shown in , such power-law modified gravity theories can serve as models for dark energy. 

Before we close this section, we need to stress that the absence of an analytic solution in the case of a singular 
bounce is exactly due to the existence of the singularity. In other bouncing cosmologies, for which no singularity 
occurs, this lack of analyticity no longer persists as a problem, and the exact behavior of the F{G) gravity can be 
found, like for example in Ref. 25|. 
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IV. STABILITY OF THE F{G) GRAVITY SOLUTIONS NEAR THE BOUNCING POINT 


The FRW equations of Eq. (flOl) for the F (G) gravity constitute a dynamical system which determines the behavior 
of solutions which satisfy these equations. Stability of a solution of this dynamical system against linear perturbations 
would mean that this solution is one of the final attractors of the theory. On the contrary if a solution of the system 
is unstable against linear perturbations, then it is obvious that this solution is not a final attractor of the theory. The 
focus in this section is on the stability of the solutions we found for the F{G), near the bouncing point, against linear 
perturbations of the solutions. What it is expected is obviously that the solutions we found near the bouncing point 
are unstable, since the cosmological evolution does not stop at the bouncing point, but continues and the Universe 
starts to expand. Therefore, the local solutions we found that describe the F{G) gravity near the bouncing point, 
should be unstable against linear perturbations, since they do not describe the whole evolution, but a small part of 
it, near the Type IV singularity. In the rest of this section, we shall use some convenient variables and we shall study 
the stability of the dynamical system of Eq. m against linear perturbations. 

We adopt the techniques and notation used in [^, so we consider linear perturbations of the solution g{N) = Fl^, 
of the following form, 

g{N) ^ g{N) F 5g{N) (29) 

with the function g{N) satisfying the ERW equations (flUl) . Expressing the dynamical system of Eq. (fTUl) in terms of 
the function g{N), we get. 


288g^iN)F''iG) [ {ig'{N)f + giN)) g" {N) + Ag{N)g'{N) + AgiN)g'{N)\ (30) 

6g{N) + F{G) - 12g{N) {g'{N) + 2g{N)) F'{G) = 0 . 

The conditions that ensure the stability of the dynamical system (1501) against linear perturbations, are the following, 


(j>0, ^i>0. 

Jl Jl 

where we introduced the variable Ji which is equal to, 

Jl = 288g{NfF''{G ), 


(31) 


(32) 


while the variable J 2 stands for, 

J 2 = A32g{N)^(^{2g{N)+g'iN))F"iG)+8g{N)(^g'iN)^+g{N){Ag'{N)+g"{N)))F"{G)) , (33) 

and in addition, the parameter J 3 is equal to, 

J 3 = 6(1 + 24g(iV)( - 8g{Nf + 3g'{Nf + 6g{N){3g'{N) + g"(iV)))F"(G) (34) 

+ 245(iV)(4g(iV) + g'iN)) (^g'{Nf + giN){Ag{N) + g"{N))) F"{G)) . 

Having these at hand, and also the stability conditions, let us investigate whether the solution (1551) is stable towards 
linear perturbations. The F{G) gravity of Eq. (l23l) can be further simplified, since a satisfies the condition 1 < a < 2, 
so the most dominant term near the bouncing point (or equivalently at the small G limit as we showed in the previous 
section), is the second term in Eq. (1231) and hence the F{G) function can be approximated by, 

F{G) ~ BG-^ . (35) 


By using the latter form of the F{G) gravity, the variables Ji, J 2 and J 3 can easily be computed and we give their 
detailed form in the Appendix B. Note that the function g{N) is expressed in terms of the e-folding number N which 
is equal to N = In a, since we have set cq = 1, so for the Hubble rate of Eq. ®, the function g{N) reads. 


9iN) 




where 7 = 3 ^. So finally, the stability conditions for the F{G) gravity of Eq. (|551) read, 


(36) 


J 2 3 {f^N{2N + 7 ) + 8N^^P*-f{-l + 4N + 27 )) 
Jl " 2/2^2 


(37) 
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and moreover J3/J1 is equal to, 


.31+1^41^iV-2+^(l - 3af/3‘^(jf 




Ji Bfoa{-l + 2 a) * ^ m/ y j2 

2 ^^ 3 ^i^-B/oA^-'>'a(-l + 2 a) (8A^2 _ 3(2 _ 3^)^) 

(l-3a)2/32(27V + 7)2 

2 ^^ 3 ^T^Ba(-l + 2 a) ( 7 )( 7 (-l + 27))^ 

iV(l-3a)2(7)2 j ■ 


(' + 


(38) 


Since we are interested in the behavior of J 2 /J 1 and J 3 /J 1 near the bouncing point, we should what the limit t —>■ ts 
means in terms of the e-folding number N. Actually, since when t ^ tg, the scale factor tends to unity, then, in the 
limit t — 7 > ts, the e-folding number tends to zero. Therefore, we shall find how J 2 /J 1 and J 3 /J 1 behave for A^ —>■ 0 
and taking the limit of the expressions appearing in Eq. (1371) and (1381) . we obtain. 


^ = 3+^ 
Ji ^ 2N' 

Jl 


(39) 


where the parameter A is positive, and can be found in the Appendix B. As it is obvious from Eq. (I39p . the parameter 
J 3 I Jl is negative and therefore the system is unstable, as we anticipated. This means that the solution (13511 is not 
a final attractor, and therefore it can describe the system near the bouncing point, but for a small time, since the 
system continues it’s evolution, after that point. What now remains is to study the evolution of scalar perturbations 
of the cosmological evolution, with emphasis on the description near the Type IV singularity. This issue is addressed 
in the next section. 


V. EVOLUTION OF SCALAR PERTURBATIONS 


Having the qualitative behavior of the F(G) gravity that generates the singular bounce, it is worth examining 
how the scalar perturbations evolve, assuming the flat ERW background of Eq. ([9|). So we consider scalar linear 
perturbations of the flat ERW background of Eq. ( 0 ), of the form, 

ds^ = —(1 -I- '0)df^ — 2a{t)dij3(ltdx'‘ + a{t)'^ {Sij + 2(j}6ij + 2didj"f) dxMx-^ , (40) 


with ijj, (j), 7 and /3 being the smooth scalar perturbations. Eor the perturbation study we follow the approach and 
master equation given in Ref. [i^, but we specify everything for the F{G) case, which is a special case of F{R,G) 
gravity studied in [d^. Perturbations are usually analyzed in terms of gauge invariant quantities, for convenience rea¬ 
sons, therefore we shall be interested in the following gauge invariant quantity (the corresponding comoving curvature 
perturbation), the evolution of which we study in this section. 


^ = (j) 


HS^ 

i 


(41) 


with f = ^- In the F{G) gravity case, the scalar modes propagating, contain no terms, so no superluminal 
propagation occurs, and only the usual terms appear [d^. The perturbation equation that governs the scalar 
modes in F{G) gravity is the following. 


a{t)^Q{t) dt 




ity 


= 0 . 


(42) 


where we can see the above equation has the usual form for scalar perturbations, in which terms dominate in the 
evolution. It is conceivable that the speed of propagation is determined by the term which for F{G) theories 

of gravity is defined to be. 


Bi{t) = l + 


2 H 

IP 


(43) 


















Moreover, the term Q{t) appearing in Eq. (1^ . for the F{G) case is equal to, 


Q{t) 


6 | 

( >Ff 

)"g2 | 

[l + AF''{G)GH^ 

( 

l + QHF''{G)G^ 

^2 


(44) 


where the prime denotes differentiation with respect to G, while the dot as usual denotes differentiation with respect 
to the cosmic time. Note additionally that we used the fact that ^ 

It is conceivable that finding an analytic solution of Eq. (l42ll is rather difficult, so either a numerical study or an 
approximate solution is required. We shall choose the latter approach and seek for an approximate solution near the 
bouncing point. Before continue, we rewrite the differential equation (1421) as follows. 


a{t)^Q{t)^ + ^3a(t)^aQ(t) + a{t)^Q(t)^ <i> + Bi(t)Q{t)a{t)k^^ = 0 . 


(45) 


After some tedious calculations, by using the resulting F{G) gravity near the bouncing point, namely the one appearing 
in Eq. (I23L and by keeping the most dominant terms near the bouncing point, the differential equation that governs 
the evolution of perturbations near t = ts reads, 

(< - - (t - = 0 , (46) 


where the parameters 17^ {i = 1,2,4) are given in the Appendix D. Note that we omitted a term ^ 
which is subdominant compared to the term {t — The parameter 173 can also be found in the Appendix D. 

Then, by solving the differential equation (HSl) we obtain the following analytic solution which describes the evolution 
of scalar perturbations near the singular point t = tg, 

$(t) = + A 2 X=(-i+“) , (47) 


with X = t — tg and the constants jj, and ( depend on the parameters 17^ and their full detailed form appears in 
the Appendix D. Also the function J^(j/) is the Bessel function of the first kind and in addition, the parameters A^, 
i = 1 , 2 , are given below, 

Ai = ("-1 + 


Note that in the parameters Aj, f = 1, 2, appear the constants G 3 and G 4 , which are arbitrary constants of integration, 
which result after solving the differential equation of Eq. (j45l) . without any initial conditions. Therefore, the solution 
ofEq. gll) is a general solution, and the constants of integration of this general solution are the parameters G 3 and 
G 4 appearing in the parameters A^, z = 1,2 of Eq. (H51) . Since we are considering the limit x —>■ 0, we can further 
approximate the solution, by using the limit of the Bessel function for small arguments, which is, 



— a' 


1 

a 

a 


‘17^ 


^17| G3r 


a 




GiV 


— 1 + 0 
172 


— 1 + 0)174 
2174 


(48) 


+ 


o174 


174 — 0I74 174 — 0I74 174 — 0I74 


Jfiiy) 




(49) 


so the approximate evolution of the scalar perturbations dUl) near the Type IV singularity behaves as follows. 


$(t) ~ A 2 


2 -M^M 


(2-2o + <»2),, 

2 (-l + c) 


(50) 


So the resulting expression for the evolution of perturbations in the absence of matter fluids is a described by a 
power-law function of the variable x = t — tg. It is worth to further investigate the power spectrum and check whether 
it is scale invariant. Note however that we already are within an approximation and therefore we should stress that 
our results should considered only approximate and also that the full F{G) solution is needed in order to answer 
the problem in a consistent way. Nevertheless, near the bouncing point it is still interesting to find how the power 
spectrum behaves within the context of the classical F{G) theory. This may indicate how the full quantum description 
of the bounce theory will remedy any problematic issues. 
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We start from the gauge invariant variable $ given in Eq. m, which as was shown in (4^, satisfies the following 
second order perturbed action, 

Sp = I dxMt)^Qs 0$ - , (51) 

where Qs = and Q{t) appears in Eq. (l44ll . Following the standard approach in perturbation theory 

the power spectrum of curvature perturbations for the scalar field d) is, 

= (52) 

It is straightforward to show that the power spectrum is not scale-invariant, by simply looking the forms of fli, A 2 
and C from Appendix D. As it can be seen, the wavenumber k is contained only in fli, A 2 and in C implicitly via fli. 
These have the following functional dependence with respect to k, 

A2 -nr^, (53) 

and since the power spectrum depends on the combination A^^, it follows from Eq. (1531) that. 


Vr 




C'4(t - ts) 


{2-2o, + o,2)^ 2 

2 (-l + a) 


k—aH 


(54) 


However, we cannot conclude at this point if the spectrum is scale invariant or not, since the parameter which is 
the constant of integration appearing in A 2 in Eq. (1481) . depends on k, and its exact form will depend on the initial 

(2-2a + <»2)M 

conditions of the vacuum form of 4>(t). In addition, the the term {t — ts) also has a dependence on k, since 

the power spectrum is computed at k = aH, that is, at the horizon crossing. So let us calculate these in detail. 

Before we start it is worth recalling that we are working for cosmological times for which t — tg —>■ 0, and hence, 
the conformal time r, defined as dr = a~^{t)dt, is approximately equal to t, since for t — ts —>-0, the scale factor 
appearing in Eq. o, behaves as a ~ 1. Also, since a ^ 1, for t — tg —>■ 0, then we have that k ~ H at the horizon 
crossing, and this means that. 


/3(t - tg)“ ~ k , 


which by solving with respect to (t — ts), yields. 


t-ts 



(55) 


(56) 


Hence, Eq. (1561) shows how t — ts behaves as a function of the wavenumber k near the horizon crossing. We now 
proceed to find the ^-dependence of (74. To do so, we shall introduce the very frequently used in the literature [4^.l43| 
canonical variable u = Zg*!*, with Zs = Q{t)a{t), and since a{t) ~ 1 for (t — ts) —>■ 0, we have that Zg ~ Q{t), and 
therefore. 




(57) 


and recall that Q{t) is defined in Eq. (HU) . In terms of u, the action of Eq. (HD) becomes near the bounce, 

Suzz I d^dr ^ --^(Vm)^- k — 

^ 2i Za 


(58) 


where the prime indicates differentiation with respect to the conformal time, which as we demonstrated earlier is 
approximately identical to the cosmological time t, for (t — tg) —>• 0. Also the action of Eq. (1551) . is defined modulo a 
factor of a~^, which is approximately equal to one near the time instance < ~ tg. The vacuum state of the canonical 
field u is the Bunch-Davies quantum fluctuating vacuum (43| at exactly t = tg, hence u ^ ^ ^ ■ Note that imaginary 

phase will disappear when we shall compute the norm of the comoving curvature |<I>(t = tg)P, so the relevant form of 
u for the k dependence issue of \s, u ^ Owing to Eq. (EZD, we obtain that, 

1 


4>(f = tg) -- Ca 


VkQ{t) 


(59) 
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By using the fact that F{G) for {t — ts) —>■ 0, is approximated by Eq. (l24l) . we get that the function Q{t) is 
approximately equal to, 


Qit) 


2 ®+-w<»3 ^+-1^+3“ - 2a)^ (a/3^) 

(1 - 3a)4^6 


and since {t — ts) is given by Eq. (|55)) . we get, 


(60) 


Q{t{k)) Zik “ = Zik 


(61) 


where Zi stands for, 

^ 2“^+^=Tfe3-l+:rrfe^2(i _ 2 q,)2 (q,/33)^T+3^ 

^ (1 - 3a)4/36;31/a ' 

By substituting Q(t) from Eq. (ICT) in Eq. (I59p . we finally obtain that that C 4 behaves as. 


(62) 


C 4 


1 1 
Z\ \fkk~'^ 


. 


(63) 


So combining Eqs. (1551) . (1551) and (1551) . we finally get that the fc-dependence of the power spectrum Vr is of the form, 


7 (2-2<. + .3)^ 

Vr - / c 2+3+ 2(- i +.) 


(64) 


and hence we conclude that the power spectrum is not scale-invariant. We can use the approximate value for the 
power spectrum of primordial curvature perturbations given in Eq. (1551) . in order to calculate the spectral index of 
primordial curvature perturbations Us as follows: Combining the expression for /r given in Eq. (11031) . and also the 
values of fl 2 and 114 given in Eq. (I102|) . we get that fi becomes = 11/(1 — a). Therefore, the spectral index of 
primordial curvature perturbations is equal to. 


^ _ d\nVn 7 , , 2-2a + a2 ^ 11 

~ d\nk 2 “^ 2(a-l) ^ 2(a-l)2- 


(65) 


Therefore, the spectral index can be in agreement with the latest observational constraints on Ug 0) for two values 
of a: one which satisfies a < — 1 and one which satisfies a > 1. The a < —1 case leads to a Big Rip singularity, while 
the a > 1 case leads to a Type IV singularity, which is the case we studied in this paper. However, the value of a 
would need to be quite large, and so it would be a large deformation of a symmetric bounce which we assumed in this 
paper, therefore this result does not appear to be very physical. 

But still we need to interpret correctly this result, since what we have at hand is that the classical approximation 
of the F{G) theory that describes the singular bounce ([5]) near the bouncing point, fails to produce a scale invariant 
spectrum. This could be a strong motivation to use a LQC corrected F{G) gravity theory and investigate 
whether the same picture persists with quantum corrections added. If the answer lies in the affirmative, then this 
could probably mean that this effect is a feature of the Type IV singularity, that needs to be further investigated. 

However, if we use more physical values of a in the range I < a < 2, this would yield a very red spectrum, and one 
can strongly doubt that quantum corrections, e.g. from LQC corrected F{G) gravity, could adjust this to a nearly 
scale-invariant spectrum. This argument is very rigid and solid and therefore should be investigated to explicitly 
demonstrate its validity, a task that we defer for a future work. 


A. Discussion 

Before we close this section, we need to discuss in detail the physical results we obtained in this section. As we 
demonstrated, the power spectrum of primordial curvature perturbations is not scale invariant, when this is evaluated 
for cosmic times near the bouncing point at < = tg, which is also the point at which the Type IV singularity occurs. 
However, we need to discuss the physical implications of this result. To this end, let us recall some fundamental 
issues for the dynamics of perturbations and related issues, but in the context of bouncing cosmology. Recall that 
the problem of initial conditions in the standard Big Bang cosmology was due to the fact that the Universe appears 
to be nearly flat and homogeneous in large scales, which cannot have causally communicated in the past. Therefore, 
a successful bouncing cosmology should in some way provide an elegant solution to these problems. 
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Before we see what happens in the bouncing cosmology we studied in this paper, it is worth to present what 
happens in inflationary theories, so that we can compare the inflationary picture to the bouncing cosmology picture. 
In the inflationary cosmology picture, the primordial curvature perturbations ^ which are of interest at present time, 
during inflation where at subhorizon scales, with horizon referring, usually in the existing related literature, to the 
Hubble radius l/a(t)H(t). Then these perturbations freeze once they exit the horizon, which happens when the 
contracting horizon becomes comparable to their wavelength. Subsequently, these freezed perturbations become 
classical superhorizon perturbations which re-enter the horizon as the horizon expands again during the Hubble 
evolution of the Universe, after the reheating. Eventually, the gravitational collapse of the freezed superhorizon 
perturbations leads to the formation of the large-scale structure of the Universe and the Cosmic Microwave Background 
anisotropies correspond to superhorizon modes which were initially subhorizon during inflation, and freezed after the 
horizon crossing. Moreover, anisotropies are due to modes which have re-entered the horizon, but still, these modes 
have the same origin on subhorizon scales during inflation. 

At this point a more concise presentation is needed so that we gain deeper insights of the full picture. The novel 
feature of the inflationary description, is that during inflation, the Hubble radius, to which we refer as horizon, actually 
shrinks dramatically. So at the initial singularity the horizon was very large, and the primordial perturbations where 
actually as subhorizon scales, since their comoving wavenumber was at subhorizon scales, that is, 

k ^ H{t)a{t). (66) 

Note that in principle perturbations are created at all length scales, but the most relevant are those which their 
wavenumber is at subhorizon scales. Also note that we switched our description by using the wavenumber, but it is 
equivalent in using the wavelength, in which case Eq. (1551) would be A <C {H{t)a{t)) During the inflationary era, 
the Hubble radius (the horizon) shrinks, so at some point the cosmologically relevant perturbations of wavenumber k 
satisfying Eq. (1661) exit the horizon and freeze, meaning that these become classical. So these become superhorizon 
perturbations, in which case, their wavenumber satisfies, 


k <C a{t)H{t). 


(67) 


Once the horizon crossing occurs, the comoving curvature perturbations corresponding to the wavenumber k cease to 
be of quantum nature, and the corresponding quantum expectation value of the comoving curvature perturbation is 
practically the classical ensemble stochastic average of a classical stochastic field. This is what we meant by freezing of 
these modes. The conservation of the average value of the comoving curvature perturbation in superhorizon scales is 
what actually enables us to relate the predictions corresponding to the time that the horizon is crossed, which actually 
corresponds to high energies, to the observable quantities corresponding to the horizon re-entry of the modes after 
reheating, which in turn corresponds to low energies. Note that the era between the horizon exit and re-entry is an era 
which is not fully understood, even up to date. Also this issue may appear also in bouncing cosmology, since in some 
cases what is needed is a quantum bounce description followed by some other model, see for example [3^ . We shall 
discuss this issue in more detail and more concretely later on in this section. Coming back to the inflationary picture, 
at present time we are able to compute the inflationary observable quantities, because the subhorizon wavelengths 
freezed out at horizon exit and evolved in a classical way until nearly the present time era, after reheating and after the 
horizon re-entering. Note that, in principle, the primordial perturbations may freeze out well before the inflationary era 
ends, so these correspond to an expanding quantum era. Hence the primordial curvature perturbations corresponding 
to the expanding quantum era may be directly related to the Cosmic Microwave Background observables and also to 
other present time observables, since the quantum fluctuations freeze after the horizon exit. 

Let us now turn our focus to the singular bouncing cosmology of Eq. ([5]). In this case, the Hubble radius as a 
function of the cosmological time is equal to. 


rnit) 




( 68 ) 


where for notational simplicity we denoted the Hubble radius as obvious from Eq. (I68L 

in the case of a Type IV singularity since a > 1, the Hubble radius at the bouncing point, which is also chosen to 
be the singularity point, the Hubble radius diverges due to the existence of the term ^ {t — Hence all the 

cosmologically relevant modes are in subhorizon scales, since k ^ H{ts)a(ts) = 0 at that point. Equivalently, at 
the singularity point, the cosmologically relevant modes are inside the Hubble radius which is infinite at that point, 


^ We are referring to the comoving curvature of Eq. gTJ, SO these are fluctuating vacuum scalar perturbations 
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so the wavelength of these modes satisfies A <C Immediately after the bouncing point, the Hubble radius drops 
and starts to shrink as the time evolves. This can also be seen in Fig. [T] where we can see that the Hubble radius 
drops after the bouncing point, in a radical way. At the left plot, we plotted until t ~ 10“^°sec and at the right 
plot until t ~ 10“^°sec, and as it can be seen, the Hubble radius rn fraction corresponding to the two cases is of the 
order — 10^^. Subsequently, as the Hubble radius shrinks, the cosmologically relevant modes will eventually 




FIG. 1: The Hubble radius th = 


order 


as a function of time, for ta = 10 


a{t)H{t) ^ wi. Awi- US J.W sec, o 13/11, j3 

be seen from the left and right plot, the Hubble radius rn fraction of corresponding to t ~ 10“^'^ 
~ 10 ^ 


ra(lO-lO) 


= 0.001(sec)“““ 
sec and t ~ 10“^'^ 


As it can 
sec is of the 


exit the horizon, when the Hubble radius becomes of the order of their wavelength X ^ ru- But which modes can 
be cosmologically relevant in the singular bounce at hand? In principle, the time era near the bouncing point is 
governed by the quantum theory of gravity, so after the bouncing point we may still have the quantum era primordial 
modes. But to which cosmological times does the “near the bounce” expression refers to? Since we assume that 
tg ~ 10“^^sec, then near the bouncing point from a mathematical point of view, corresponds to cosmological times 
of the order t ~ 10“^°sec, which can also be considered as being near the bounce, since what we assumed in the 
calculation of the spectrum of the primordial curvature perturbations is that t — 0. Thereby, in inflationary 

terms these cosmic times correspond to times after the exit from the inflationary era, so the relevant for today modes 
have already exited the horizon well before t ~ 10“^°sec. Hence for the singular bounce, the modes we took into 
account in the calculation of the spectrum of the primordial curvature perturbations, are actually the cosmologically 
relevant for today, quantum modes. Therefore, in the singular bounce case, the modes with wavelengths of the order 
of the horizon corresponding to cosmic times near the bouncing point, are the cosmologically relevant, since these 
can reveal the quantum era of primordial expansion. After these modes exit the Hubble radius, freeze out and the 
quantum expectation value of the comoving curvature perturbation is described by the classical ensemble stochastic 
average of a classical stochastic field. Practically, the conservation of the average value of the comoving curvature 
perturbation at superhorizon scales will eventually enable us to relate the horizon crossing predictions (high energy 
ones) to the late-time ones, which correspond to the horizon re-entry era (low energy). 

However we need to explicitly demonstrate that in the case of the singular bounce we studied in this paper, the 
comoving curvature perturbations are conserved after the modes exit the horizon. This is important, since it is not 
granted that the comoving curvature perturbations will be conserved, like for example in the inflationary cosmology. 
An example for which the comoving curvature perturbations grow after the modes exit the horizon is the matter 
bounce cosmology case (d^j . so now we study the evolution of the comoving curvature perturbations in the context 
of the singular bounce. We are interested in cosmological times for which the Hubble rate and scale factor satisfy 
k <C a{t)H{t), so for times much more later than the horizon crossing. Note that this does not mean that t corresponds 
to late-times, or t ^ 1, but the cosmic time is t ^ with tn the time at which the horizon crossing occurs. Since 
k <C this means that the last term in the differential equation of Eq. (H^ . can be neglected, so the differential 

equation becomes, 

(69) 

which can be easily solved and the solution is, 

$(t)=Ci+C2/^(^)3 ^(^)dt, (70) 

with Q{t) being defined in Eq. (j44p . So in order to see if the comoving curvature perturbations are conserved after 
the horizon crossing, we need to examine the behavior of the integral term in Eq. GOl). Obviously, the key point to 
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determine the behavior of the comoving curvature perturbations, is to determine Q(t), and therefore to find which 
F{G) gravity describes the cosmological evolution at the cosmological times for which k ^ a{t)H{t). In order to 
proceed we need to classify the problem in two subcases or scenarios more preferably, which we list below: 

1. In the first scenario, which we call Scenario I, the cosmic times for which the relation k <C a(t)H{t), satisfy 
t ^ tg, and also t <C 1, where tg the time at which the Type IV singularity occurs. This is the most plausible 
scenario, from a physical point of view, and it could be realized like this: Suppose that tg = 10“^®sec and the 
time long after the horizon crossing is at t = 10“^®sec. Hence, in this case, the time t is 10^° times larger than 
tg and also satisfies t <C 1. Note that the time at which horizon crossing occurs, is somewhere in between tg and 
t, that is, tg < tH < t. 

2. In the second scenario, which we call Scenario II, we have again t ^ tg, but t > 1. 

The Scenario I is much more likely to occur, since the time for which k <C a{t)H{t) holds true, is after the horizon 
crossing, which for inflationary cosmology is (possibly) of the order 10“^°sec, so this is the most appealing case and 
we start our analysis with Scenario I, but for completeness we also deal with Scenario II later on. 

So in the context of Scenario I, the cosmic time satisfies t tg, but still t <C 1, and consequently, in this case 
t — tg —>■ 0. Therefore, the F{G) gravity that is responsible for the cosmological evolution at the time t, is still given 
by Eq. (1^^ . and therefore, the function Q{t) can easily be calculated to yield, 

Q{t) ~ , (71) 


where in Eq. dn]) we kept the most dominant term, and also we used the fact that t — tg ~ t <C 1. Therefore, since 
a{t) ~ for the singular bounce (see Eq. ([T|)), the term behaves as follows. 


1 


f 


4a 


a(t)^Q(t) 

Clearly, the exponential dominates, so the integral term decays as t increases, that is, 

/ a{t)3Q{t) 

and hence, the comoving curvature perturbation is approximately equal to. 


(72) 


(73) 


$(t) = Ci, 


(74) 


which means that the comoving curvature perturbation for the Scenario I is conserved after the horizon crossing. 

Now we turn our focus on Scenario II, for which t ^ tg and t > 1. As we already mentioned, this scenario probably 
corresponds to an era much more later than the bounce, and possibly much more later than the horizon crossing, so 
it is rather less physically appealing. Regardless, we shall study this scenario for completeness. It is conceivable, that 
the F{G) gravity which describes the cosmological evolution is no longer given by the one appearing in Eq. (IMl) . since 
{t — tg) is not small anymore. Therefore we need to find the F{G) gravity that describes the cosmological evolution. 
This is however a formidable task, since by employing the reconstruction method we used in the previous sections, 
we end up to the following differential equation, 

2t“+ip(t) - 2f+^°‘al3‘^P{t) + 1 = 0, (75) 

which is very difficult to solve analytically. So in order to proceed, we speculate about the possible behavior of the 
F{G) gravity. If for example the F{G) gravity has polynomial form, that is F{G) ^ B G^, then, since for large t, the 
Gauss-Bonnet invariant becomes approximately equal to G ^ the resulting expression for Q{t) is, 

Q{t) ~ 2 : 3 ^ 2 +®“+®“^', (76) 


with Z-i being equal to. 


Z3 = 212+6732-^2752^2^8-^87 ^ 


(77) 


Therefore, the term ^ l/(a(t)®(5(t)), is equal to, 

1 I 


a(t)3Q(t) e/o*“+'t-2+8a-S8a7 ' 


( 78 ) 
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Hence, even in the case that 7 is a large negative real number, the exponential in the expression appearing in Eq. 
(1751) dominates and thereby the integral in Eq. (1701) decays and becomes subdominant. Hence in this case too, the 
comoving curvature perturbations after the horizon crossing are conserved, since $(t) = Ci. 

In the case that the F{G) is such so that Q{t) dominates over the exponential scale factor, then the integral in Eq. 
(1701) dominates the evolution, and possibly the comoving curvature perturbations might grow as the time passes. For 
example, if F{G) ^ e~'^ , then the function Q{t) becomes approximately equal to, 

(79) 

with Ai = and therefore the term ^ 1 /{a{t)^Q(t)) becomes approximately equal to, 

a(t)3Q(t) ^ 

which clearly does not decay. So the integral in the expression (|70l) dominates the evolution of the comoving curvature 
perturbations after the horizon crossing, and the perturbations grow as the time grows larger. Therefore, in the case of 
the Scenario H, the curvature perturbations depend strongly on the form of the F(G) gravity. However, the Scenario 
H is rather unlikely to occur, since the requirement t > 1 sec means that the Universe is at the lepton epoch, which is 
much more later than the time that the horizon crossing occurred. At the time t > 1, it is possible that the singular 
bounce does not describe the Universe anymore, since another description must be found to generate the Hubble 
radius expansion, since the Hubble radius decreases in the context of the singular bounce. It is highly likely that 
a scenario like the one used in [s^, will take place and describe the Universe’s evolution. Work is in progress for 
realizing such a scenario, but it is worth analyzing this a bit more. 

The problem with the singular bouncing cosmology occurs, because there exists no mechanism to make the Hubble 
radius increase eventually, unless we assume that the singular bounce governs early-times and after some time instance, 
the evolution is described by another scale factor, as in for example in Ref. [s^, see also [s^ for an F{R) description 
of the model of Ref. [13 ■ In the case of the model studied in Ref. [131, the quantum radiation era was followed 
by a perfect fluid evolution, see [13|- So in the singular bouncing case, it is necessary to find another scenario that 
will describe the evolution of the Universe after the quantum bouncing era, for which new scenario, the re-entering 
of the modes in the expanding horizon will be possible. However, the in-between the two horizons era, leaves a gap 
in our description, and as we noted, this also occurs in the standard inflationary cosmology. This task exceeds the 
purposes of this paper, but we hope in a future work to provide a model with two or more phases, for which a 
successful cosmological description may be achieved. However, this will not be in the context of F(G) gravity, since 
as we demonstrated the spectrum of primordial curvature perturbations is not scale invariant, and this is a rather 
discouraging feature, but perhaps the LQC-corrected F{G) gravity may have more appealing features. 


VI. SINGULAR BOUNCE FROM MIMETIC F{G) GRAVITY 

As a final study, we shall investigate which F{G) gravity can generate the singular bounce of Eq. ([5]), but in the 
context of mimetic F{G) gravity. For a detailed account on this issue, see [s^l- The mimetic F{G) gravity a ppr oach 
uses the same action as the one appearing in Eq. ( 0 ), but the Jordan frame metric is parameterized as follows (47H^ . 

= -g^'^dp(j)da(l). (81) 


Upon varying the metric tensor, we get, 

= g'’"'5grug"'''dp(j)da(t>gpy - g'’'^dp(j)da(j)8gp,y- 


-2gP'^dp(j)d„5(l)gpi,. 

and upon variation of the Jordan frame action ([5]), with respect to the redefined metric g^^i and with respect to the 
mimetic scalar we obtain the following equations of motion, 




( 82 ) 


-b8 


R 

R^pya RpyQap RpaQyp, Rp,y9ap RpaQyp ^9p,y9<jp Fg (^FqG F(^G')') Qp,y~\~ 
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^d^(t>d,ct> (^-R + 8 (^-Rp, + ]^Rgpa^ V^V^Fg + 4(FgG - F(G))^ = 

= Tpi, + dp4>d„4>T, 

with Fq standing for Fq = dF{G)/dG. Moreover, upon variation of the action ([S|) with respect to the mimetic scalar 
field (j), we get, 

(^-R + 8 (^-Rp, + ^Rgpa^ V^V^Fg + 4(FgG - F{G)^ - T)^ = 0. (83) 

Since the following relation holds true [s^, 


g'"''dp(j)d^(l> = -1 


(84) 


and owing to the fact that the mimetic scalar (f> depends only on the cosmic time, we get the constraint (j) = t. Hence 
the {t, t) component of the expression appearing in Eq. (l82l) becomes. 


2H + + 1&H(H + + 8H' 

at 


r 2 d^FG 


dF 


- {FgG - F{G)) = -p. 


Upon integration of Eq. (1551) . we get. 


-R + 8 



F^Rg 



S/PV^Fg + 4(FgG - F(G)) + p - 3p 



(85) 


( 86 ) 


which can be rewritten as follows. 


H + 2H^ + AH 


2 , .TT2d'^FG 


dF 


AH 


(2H + 3i7^) 


‘"^=G(FoG-F(G)) + U? = - 


where G being an arbitrary constant. 

The combination of Eqs. (IMI) and (1^ . results in, 

pd^Fg , ■ „,dFg 1, , G 

H + + 4H(2H - 

For convenience, we introduce the function g{t), which is defined to be, 

dFG 


git) = 


dt 


G 


and satisfies the following equation, 

4 ^ 2 ^ + ^H{2H - H^)g{t) = -H-\{p + p)-^. 
at 2 a'^ 

Then, the general form of the solution of the differential equation ([89|) is, 


git) = flo 




4iJ2 


dtiB{ti) exp 



(87) 


( 88 ) 


(89) 


(90) 


where Bit) = —H — \ ipF p) — ^ and in addition po is an integration constant, while, Hq = i7(0). Having Eq. (ITOl) . 
for a given cosmological evolution in terms of the Hubble rate, we can easily obtain the function git). Consequently, 
the function F[t) reads. 


FGit) = j git)dt. 

Then by exploiting the expression ([8]), and solving with respect to t, we have the function t = t{G). By substituting 
this to FG{t), and by integrating with respect to G, we easily obtain the F{G) function. Let us apply this method for 
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the case of the singular bounce of Eq. ©■ By substituting the Hubble rate of Eq. in Eq. O, and by keeping 
the most dominant terms, we obtain, 


git) 






4/2(1+ a )2 



(91) 


Integrating this with respect to t, we obtain the approximate form of the function Fait) near the bouncing point, 
which is. 


Fait) = 


{t - ts)^ {-16goF[Q (3 - 5q; + 2a‘^) + + t{ts - 2tsa) + (l - So; + 2a'^)) (2C + p + p)) 


16/2(-1 + a)(l + a)2(-3 + 2a){-l + 2a) 


(92) 


/o 


l+a 


(l + a)(2-5a + 2a2) ' 

Using ®, and also Eq. dini), we obtain the function Fg{G), so upon integration with respect to G we finally get. 


F{G) ~ -24/3a(l + a)3 (a7G"T^) 
-24/((a(l + a)3 oi ^ 


(93) 


06 


G-l + 3a 02 + 3 + 04 ts + 24 i-3“G-1+3“ - 


+ a)3 


-24/3a(l + a)3(2G + p + p)o5aiG-i+3- ( 


06 


ts + 241-30 G-i+3“ - 


/3a(l -k a)3 


where the parameters o^, j= 1,2,...7 are given in Appendix C. Since we are considering the limit t ^ tg, which means 
G —>■ 0 as we explained earlier, by keeping the most dominant term in Eq. (IM)l . we get. 


F(G) 


24/ga;(l + q) 3(2G +p + p)t2a5ai ^^^ 
06 


(94) 


It is conceivable that the resulting mimetic F{G) gravity of Eq. (IMl) is different for the vacuum F{G) gravity of Eq. 
(|2^ . but we need to further analyze this issue. In the context of the mimetic F{G) gravity, some extra conformal 
degrees of freedom arise in the FRW equations of motion. Therefore, in the context of mimetic F{G) gravity, we have 
a new reconstruction method in which we can choose the internal degrees of freedom and a specific F{G) gravity so 
that some fixed cosmological evolution is generated. This is different from the ordinary vacuum F{G) gravity case, 
since in this case no internal degrees of freedom are taken into account, so the F(G) gravity that can generate the 
same Hubble rate as the mimetic F{G) does, is in principle different from the resulting expression of the mimetic 
F(G) gravity. Of course in both cases we are using approximations, so one should be cautious when dealing with 
both theories. Finally, let us note that in the mimetic F(G ) gr avity, much more freedom is offered for successfully 
generating various cosmological scenarios, see for example [50j. This is because of the presence of these internal 
conformal degrees of freedom. It is then easy to reconstruct any cosmology by suitably adjusting these degrees of 
freedom, and in some approaches the potential and the Lagrange multiplier. It is questionable however if these results 
can be trusted, because the resulting picture is complicated. If simplicity is to be an important feature of a physical 
theory, then probably these theories are of mathematical importance only. However, concordance with observations is 
an appealing feature, and therefore these theories can be valuable from a physical point of view. Since this discussion 
should be addressed in more detail, we defer it to a future work. Regardless of the differences, in both cases, the F{G) 
gravity that can realize the singular bounce of Eq. ®, is described by a power-law function. But in our opinion, the 
vacuum F{G) is conceptually a simpler theory, so from this aspect it offers a more appealing physical description of 
the singular bounce. 


VII. CONCLUSIONS 

In this paper we studied a bounce cosmology with a Type IV singularity at the bouncing point, in the context 
of classical F{G) gravity. Particularly, we investigated which classical pure (vacuum) F{G) gravity can generate 
the Type IV singular bounce cosmology, emphasizing for cosmic times near the bouncing point. As we explicitly 
demonstrated, the resulting F{G) gravity has the form F{G) ^ G 2 G + B G , so it is a power-law modified gravity 
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theory. Since this result holds true only near the singularity point, we discussed the possibility that this F{G) gravity 
is the limiting case of some viable F{G) gravity, in which case the full solution would also be interesting, since the 
late-time and early-time acceleration could be simultaneously described by the same theory. We also discussed the 
stability of the resulting F{G) theory, from a dynamical point of view, examining if it can be the final attractor of 
the theory. As we anticipated, the answer to this question does not lie in the affirmative, and hence instability of the 
solution cannot be avoided. This feature is welcome, since the cosmological evolution does not stop at the bouncing 
point, and therefore the resulting F{G) gravity was not anticipated to be a stable solution of the cosmological system. 
Moreover, we investigated how the scalar cosmological perturbations of the background flat FRW metric, behave near 
the bouncing point, and we explicitly calculated the spectrum of primordial curvature perturbations. As we showed, 
the spectrum is not scale invariant and as we claimed in the main text, this result should be further investigated. 
This is due to the fact that we cannot be sure if it is a universal feature of the theory that owes its existence in the 
Type IV singularity, or it is an artifact of the approximations we made to obtain the resulting F{G) gravity. The 
latter seems more plausible, however this feature has to be thoroughly addressed. Another important point that we 
need to stress, with regards to the non-scale invariance of the spectrum of primordial curvature perturbations, is that 
since we are studying a classical theory near the bouncing point, it might be possible that at these cosmic time scales, 
quantum effects take place. So effectively the lack of scale invariance in the power spectrum might be an effect of 
the classical ap proach to the problem, so the same problem should be addressed in the context of Loop Quantum 
Cosmology [sy, and particularly in the context of F{G) LQC, which was developed in [i^. We also studied which 
mimetic F{G) gravity can describe the singular bounce near the bouncing point, by adopting the formalism of [s^ 
and the resulting F(G) gravity has a power-law functional form. 

Finally, with regards to the classical A(G) gravity approach, since the F{G) gravity is a special case of the most 
general class of F{R,G) theory Elm, [ 53 . ISSj . the same problem we investigated in this problem should be 
addressed in the context of F{R, G) theory. Actually, this problem should also be compared with the F{R, G) gravity 
inflation properties, as was done in [s^, but this time by using a Type IV singular bounce. In addition, a compelling 
task is to include matter fluids in the theory and investigate how the physical picture is affected by the presence of 
matter. Moreover, as was demonstrated in Ref (35j | that a Type IV may play a crucial role in the graceful exit from 
inflation, but the study was focused on scalar field models. It is worth to examine the effect of the Type IV singularity 
on Jordan frame F{G) theories but also in F{R) and F{T) theories. For a related work with regards to the latter, 
see [S^. We hope to materialize these projects in a future work. 
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Appendix A: Analytic Form of Q(t) and of A, B 


In this Appendix we quote the exact form of the function Q{t) and of A^B appearing in Eqs. (1171) and (1211) . 
Particularly, the function Q{t) reads. 


Q{t) = -2{t - - 72{t - 

/ {t - ((t - A)-^+“a - 2{t - 

V 2(-l-f a)(-l-f 2a)/3 

(t - G)-2“(1 - 2a) ((t - tsT - 2{t - tsYa - 2Gi^ + 2Gia/3) 


24(t-G)3“/33(. 


2{-l + a){-l + 2a)/3 
{t - - 2a) {{t - A)“^+“a - 2{t - A)“^+“a2) 
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2( —1 + o;)(—1 + 2(y.)/3 

{t — — 2a)a {{t — ts)^ — 2{t — ts)^a — 2Ci(3 + 2Cia(3) 

(—1 + oi)(—1 + 2o.)j3 


)■ 


(95) 
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Also, the exact analytic form of the coefficients A and -B appearing in Eq. (ED is the following, 
A= 11 24^/3^ (a/3^)'^ 

1_2„ 1 2„ 1 I 3 ( 1 - 20 .) l-2„ 1 - 2 ^ \ l-2e. 

. 1-42- ^ r. ^ 442. ^ 2 i-3.» 31 - 3 ° , 241 - 3 ° a \ 


B = 


24 1 - 3 ° Cl - 241 - 3 ° Cia- 


(a/3^)' 


1 - 3a + 2a2 

_ 24^+T^CiG^^"Tf3^r/33 


(96) 


Appendix B: Exact Form of the Parameters Ji, J 2 , J 3 and of A 
Here we quote the exact form of the parameters Ji, J 2 , J 3 and of A. Particularly, the parameter Ji is, 


2l+^Tq^3irT+3SBiVT'a^/3 


2 fl 2 


/o 


Jl — 


/o(-l + 3a)2(M + iv^y 


21 +^Tf 3 S 3 
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~?5~+ to ) \ 
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Also the parameter A which appears in Eq. dMD, is equal to. 
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Appendix C: Detailed Form of the Parameters ai, a2, as, 04, 05, ae, 07 


The detailed form of the parameters a^, i = 1,2, ...7 that appear in Eq. (IMI) are. 


ai 

02 

03 

04 

05 

06 

07 


( /oMl + «)0 ) ’ 

-l6goH^ (18 - 15a - lOa^ + + 2a^) 

tl (11 — 6 a + 3a^) , 
t, (5 - 6 a - Qa^ + 2a^) 
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16/o^(l + a)3(-l + 2a) (-2 + a + a^) (-9 + 3a + 2a^) 

fo^goH^o ( 24 ^ 

6 a — 9a^ — 9a^ + 6a^ 



( 101 ) 


Appendix D: The Parameters fli, Sis, Us, fl4, p and ( 


Here we quote the detailed form of the parameters Hi, H 2 , H 3 , H 4 and Ai, i = 1,2, /r and C- Particularly, the 
parameters Hi, H 2 , H 3 , H 4 appearing in Eq. (l46ll are equal to. 


Hi = 
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( 102 ) 


H 3 = 
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In addition, the parameters g and C, appearing in Eq. 621 ) are equal to, 

H 2 + H 4 
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